Non-perturbative renormalization-group approach to the Bose-Hubbard model 
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We present a non-perturbative renormalization-group approach to the Bose-Hubbard model. By 
taking as initial condition of the RG flow the (local) limit of decoupled sites, we take into account 
both local and long-distance fluctuations in a nontrivial way. This approach yields a phase diagram 
in very good quantitative agreement with the quantum Monte Carlo results and reproduces the two 
universality classes of the superfluid-Mott-insulator transition with a good estimate of the critical 
exponents. Furthermore, it reveals the crucial role of the "Ginzburg length" as a crossover length 
between a weakly- and a strongly-correlated superfluid phase. 
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Introduction. In the last two decades, the non- 
perturbative renormalization-group (NPRG) approach 
has been successfully applied to many areas of physics, 
from high-energy physics to statistical and condensed- 
matter physics (for reviews, see Refs. [ , ]). Although 
the RG is often seen as a powerful tool to study the low- 
energy long-distance physics in the framework of effective 
field theories, it has recently been shown that the NPRG 
also applies to lattice models and enables to compute 
not only universal quantities (critical exponents) but also 
non-universal quantities (such as phase diagrams, transi- 
tion temperatures, order parameters) which strongly de- 
pend on the microscopic parameters of the model (lattice 
type, strength of the interactions, etc.). This implemen- 
tation of the NPRG is referred to as the lattice NPRG [ ]. 

In this paper, we report a NPRG study of the Bose- 
Hubbard model [ ] . This approach yields a description of 
the superfluid-Mott-insulator transition which takes into 
account both local fluctuations (which drive the Mott 
transition and determine the phase diagram) and critical 
fluctuations in a nontrivial way. By comparing with the 
numerically exact lattice quantum Monte Carlo simula- 
tion (QMC), we show that the NPRG yields remarkably 
accurate results for the phase diagram. Moreover, con- 
trary to the QMC simulation, we obtain the critical be- 
havior at the Mott transition and recover the existence of 
two universality classes [ ]. We also emphasize the cru- 
cial role of the Ginzburg length as a crossover length 
between a weakly- and a strongly-correlated superfluid 
phase. 

The non-perturbative RG. The d-dimensional Bose- 
Hubbard model is defined by the (Euclidean) action 
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where tpri^) is a complex field and r G [0, /3] an imaginary 
time with f3 ^ oo the inverse temperature, {r} denotes 
the N sites of the lattice and (r, r') nearest-neighbor sites. 



U is the on-site repulsion, t the hopping amplitude and /i 
the chemical potential. (We take h = = ^ throughout 
the paper.) 

The strategy of the NPRG is to build a family of mod- 
els with action Sk = S -\- ASk indexed by a momentum 
scale k varying from a microscopic scale A down to 0. 
This is achieved by adding to the action (1) the term 

ASk = /o^<^^I]qV^q^/c(q)V^q (V^q IS the Fouricr trans- 
form of ipr)^ where 

Rk{q) = -ZA,ktk^sgn{t^){l - ^q)6(l - ^q), (2) 

with tq = — 2t ^^^-^ cos^i, ?/q = {2dt — |tq|)/t/c^ and 
Q{x) the step function (we take the lattice spacing as the 
unit length). The /c-dependent constant ZA,k is defined 
below. Since i?/c=o(q) = 0, the action Sk=o coincides 
with the action (1). On the other hand, for A: = A = 
V2d^ Ra{^) = (we use Za^a = 1) and Sa = S -\- 
ASa corresponds to the local limit of decoupled sites 
(vanishing hopping amplitude), a limit which is exactly 
solvable. For small /c, the function i^/c(q) gives a mass 
~ k'^ to the low-energy modes |q| < k and acts as an 
infrared regulator. 

The Bose-Hubbard model (with action Sk=o) can be 
related to the reference model (with action Sa) by a RG 
equation. We consider the scale-dependent effective ac- 
tion 

Tk[rA] = -lnZ,[J*,J]+ / dr^(j;(/)r + c.c.) 

[(/)*,(/)], (3) 

defined as a (slightly modified) Legendre transform which 
includes the explicit subtraction of ASk [^* , ^] • Here 
Z/e[J*, J] is the partition function, Jr(r) a complex ex- 
ternal source which couples linearly to the bosonic field 
and (j)r{T) = ^lnZ/c[J*, J]/SJ^{r) is the superfluid order 
parameter. The variation of the effective action with k 
is governed by Wetterich's equation [5], 
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where ^ is the second-order functional derivative of 
F/e. In Fourier space, the trace in (4) involves a sum over 
momenta and frequencies as well as the two components 
of the complex field (j). The initial condition of the RG 
equation is 



Fa [(/)*, 0] = Fioc[0*,0] 



Jo 



r(q)tq0(q), (5) 



where Ti,,[4>* ,<j>] = -lnZ^[J*,J] + drY^^J^'^r + 
c.c.) is the Legendre transform of the free energy 
— In Za [J* , J] of the reference system corresponding to 
the local limit of decoupled sites. The effective action Fa 
reproduces the strong-coupling RPA theory of the Bose- 
Hubbard model, which treats exactly the on-site repul- 
sion but takes into account the inter-site hopping term 
in a mean-field type approximation [6-8]. The strong- 
coupling RPA theory describes qualitatively the phase 
diagram but is not quantitatively accurate and breaks 
down in the critical regime near the superfiuid-Mott- 
insulator transition. In the NPRG technique, fluctua- 
tions beyond the RPA are included by solving the flow 
equation (4). Since the starting action ^a is purely lo- 
cal, our approach is to some extent reminiscent of various 
t/U expansions of the Bose-Hubbard model [ ]. 

We are primarily interested in two quantities. The 
first one is the effective potential defined by Vk{n) = 
^F/c[0*,(/)] with (j) a constant (i.e. uniform and time- 
independent) field and n = Its minimum determines 
the condensate density no,/c and the thermodynamic po- 
tential (per site) V^^k = ^/c(^o,/c) in the equilibrium state. 
At the initial stage of the RG, Va(^) = Moc(^) — 2dtn^ 
where Moc(^) is the thermodynamic potential in the local 
limit. 

The second quantity of interest is the two-point ver- 
tex F^ ^ which determines the single-particle propagator 

f2) — 1 

Gk = — F^ ^ and therefore the excitation spectrum. 
Because of the U(l) symmetry of the action (1), the two- 
point vertex in a constant field takes the form 

-(2) 



(6) 

in Fourier space, where q = (q, iu) and is a Matsubara 
frequency. Here (0i,02) = A/2(Re((/>), Im(^)), n = = 
^(02 +^2) ^ij is antisymmetric tensor. In order to 
solve the flow equation (4), we use a derivative expansion 



of F 
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TA,k{q; n) = VA,k{n)uj'^ + ZA,/c(^)eq + Vl{n), 
FB,fe(g;n) = F,^'(n), 
rc,fe(<7;n) = Zc,kin)uj, 



(7) 



where Cq = + 2dt (eq tq^ for |q| <C 1). This deriva- 
tive expansion is similar to the one used in continuum 

models [ - ], but the initial conditions at scale k = A 

(2) 

are here obtained from F^^"^ [Eq. (5)] and therefore al- 
ready include on-site quantum fluctuations. To reduce 
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FIG. 1. (Color online) Top panel: Phase diagram of the two- 
dimensional Bose-Hubbard model showing the first Mott lob 
with density n = 1. The NPRG result is shown by the (blue) 
points and the QMC data of Ref. [ ] by the dashed line. 
Inset: density n vs for different values of t/U. Bottom 
panel: Transition lines obtained from three different approxi- 
mations (see text). The best one (iii) is also shown in the top 
figure. 



the numerical effort, one can further approximate Va,/c(^) 
by VA,k = VA,k{no,k) (and similarly for ZA,k{n) and 
^c,/c(^)) and expand the effective potential to quadratic 
order about its minimum, 



Vk{n) 



Vb,/c + ^(n - no,/c)^ if no,/c > 0, 
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In the superffuid phase, Eqs. (7,8) yield a gapless mode 
uo = c/c|q| with velocity 



ZA,kt 



1/2 



(9) 



and a superffuid stiffness (deflned as the rigidity wrt 
a twist of the phase of the order parameter) ps^k = 
2tZA,k'^o,k- All physical quantities of interest can now be 
obtained by solving the flow equation (4) together with 
Eqs. (6,7) (and, possibly, Eq. (8)) [ ]. In the following, 
we focus on the two-dimensional Bose-Hubbard model; 
the three-dimensional model will be discussed elsewhere. 

Phase diagram. For given values of t, /7 and /i, the 
ground state can be deduced from the values of the 
condensate density no,/c=o (^o,fc=o > in the super- 
fluid phase), while the density is obtained from n = 
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FIG. 2. (Color online) NPRG flows in the weakly- correlated 
superfluid phase, t/U = 10 and kc ^ kh ^ ^ (top), and 
the strongly-correlated superfluid phase, t/U 0.062 and 
kh ^ kc ^ (bottom), for a density n — 1. The insets show 
Cfc, no,fc and ps,k vs ln(A/A;). 



— ^VQ^k=Q' It takes only a couple of seconds (depending 
on the approximation scheme, see below) to solve numer- 
ically the NPRG equations (for U and /i fixed) on a 
standard PC, so that the full determination of the phase 
diagram requires at most an hour. 

Fig. 1 shows the phase diagrams obtained from three 
different approximations: i) the effective potential Vkiji) 
is truncated [Eq. (8)] and the n dependence of Va,/c(^)7 
^A,/c(^)7 Zc,k{'^) is neglected as explained above; ii) the 
full n dependence of Zc^kij^) is included; iii) the full n de- 
pendence of Vfc(n) and Zc^ki^) is included. By including 
more functions into the analysis (i.e. going from (i) to 
(iii)) we observe a nice convergence of our results, which 
we therefore expect to be close to the exact ones, with a 
typical error, estimated from the difference between (ii) 
and (iii), roughly of order of 3%. This expectation is con- 
firmed by a direct comparison to the QMC data [ ]: the 
tip of the Mott lob {t/U = 0.060, /i/U = 0.387) differs 
from the QMC result only by (1.5%, 4%). It should be 
noted that the accuracy of the NPRG (within similar ap- 
proximation schemes) in computing non-universal quan- 
tities (phase diagrams and thermodynamics) has been 
reported in other contexts, in particular in classical spin 
models [ ] and finite-temperature field theory [15]. 

Superfluid phase. In the weak-coupling limit, we re- 
cover the results of previous NPRG studies in continuum 
models [10-12]. The strong-coupling RPA (the initial 
condition of the RG flow) is equivalent to the Bogoli- 
ubov approximation when <C 1 [ ] . The condensate 
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FIG. 3. (Color online) Top panel: Dimensionless con- 
densate density no,fc and coupling constant Xk vs ln(A//c) 
at the XY critical point (the inset shows In Zc,fc(no,fc) vs 
ln(A//c)). Bottom panel: (running) anomalous dimensions 
r]k = -/c^fc lnZA,fc(no,fc) and r]v,k = -/c^fe In VA,fe(no,fc). The 
dynamical exponent is defined by 2;^ = (2 — ?7/c + ?7y,fc)/2. 



density no,/c, the superfluid stiffness ps^k or the Goldstone 
mode velocity Ck vary weakly with k and are well approx- 
imated by their Bogoliubov estimates no, a — n, ps^A — 
2tn and ca ^ (2/7tn)^/^ (Fig. 2). On the other hand, 
from the strong variation of A^, Zc^k and Va,/c with /c, we 
can distinguish two regimes separated by the characteris- 
tic (Ginzburg) momentum scale kc = S.g^ ^ \/^{U /t)^: 
i) a (perturbative) Bogoliubov regime k ^ kc where 
\k - >^A - U, Zc,k - Zc,A - 1 and VA,k - 0. The 
spectrum crosses over from a quadratic dispersion to a 
linear sound-like dispersion at the (healing) momentum 



scale kh = — \/nU /t defined by no^k^k — ZA^ktk"^- 
ii) a (non-perturbative) Goldstone regime k <^ kc where 
A/c, Zc,k ^ k vanish with k ^ and Va,/c — takes a 
finite value. This regime is dominated by phase fluctua- 
tions, and characterized by the vanishing of the anoma- 
lous self-energy T^an,k{Q = 0) = Xk^o^k ^ k and the di- 
vergence of the longitudinal propagator G||,/c(^ = 0) = 
l/(2A/cno,/c) [16]. A weakly-correlated two-dimensional 
superfluid is characterized by kc/k^ ~ U /t <C 1; al- 
though the Bogoliubov theory breaks down at low en- 
ergy, it applies to a large part {kc ^ |q| ^ k^ <C A) of 
the spectrum where the dispersion is linear. 

As t/U decreases, the ratio kc/kh increases and even- 
tually becomes of order one (with k^ ^ kc ^ A.). In this 
strongly-correlated superfluid phase, there is no Bogoli- 
ubov regime any more and the condensate density no = 
no,/c=0 7 as well as the superfluid stiffness ps = ps,k=o^ is 
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TABLE I. Critical behavior at the superfluid-Mott-insulator 
transition and infrared behavior in the superfluid and Mott- 
insulator phases. The stared quantities indicate nonzero 
fixed-point values and 77 denotes the anomalous dimension 
at the three-dimensional XY critical point. ZA,k and VA,k 
stand for ZA,k{no^k) and VA,fc(no,/c). 



strongly suppressed (Fig. 2). 



Critical regime. Our approach recovers the two uni- 
versality classes of the superfluid-Mott-insulator transi- 
tion [4]. Away from the tip of the Mott lob, the transi- 
tion is mean-field like (with logarithmic corrections) with 
a dynamical exponent z = 2. At the tip, it belongs to 
the universality class of the three-dimensional XY model. 
Fig. 3 shows the RG flows of the dimensionless coupling 



= k^'VX'J^ {ZA^ktk 



Zc,k{^o,k) = {VA,kZA,ktk'^ 



1/2 
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1/2 



^C,/c(^0,/c), 



(10) 



when the system is at the XY critical point. The 
plateaus observed for the dimensionless condensate den- 
sity no,/c and coupling constant A^, as well as for the (run- 
ning) anomalous dimensions r]k = —kdk In ZA,k{'^o,k) and 
r]v,k = — ^Sfc In VA,/c(^o,fc), are characteristic of critical 
behavior. We find the critical exponent u = 0.699, the 
anomalous dimensions r] = 0.049, rjv = ^(1 — rj/A) = 
0.049, and the dynamical exponent z = {2 — r]-\- r]v)/2 = 
1.000, to be compared with the best known estimates 
u = 0.671 and 77 = 0.038 for the three-dimensional XY 
model [ ]. Table I summarizes the infrared behavior 
of the two-dimensional Bose-Hubbard model. Note that 
both in the superfluid phase and at the XY critical point, 
the infrared behavior is characterized by a relativistic 
symmetry {Zc,k ^ for /c ^ 0). 

Conclusion. The excellent agreement between our re- 
sults and the QMC data [1^] shows that the lattice 
NPRG, first introduced in Ref. [3] for classical systems, 
is a very efficient method to study the Bose-Hubbard 
model. This RG approach, which is implemented in mo- 
mentum space, takes into account local fluctuations and 
Mott physics while being able to describe critical fluctu- 
ations. 

We would like to thank B. Delamotte for useful discus- 
sions and B. Capogrosso-Sansone for providing us with 
the QMC data shown in Fig. 1. 
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